Universality in Three-Frequency Resonances 
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We investigate the hierarchical structure of three-frequency resonances in nonlinear dynamical systems with 
three interacting frequencies. We hypothesize an ordering of these resonances based on a generalization of the 
Farey tree organization from two frequencies to three. In experiments and numerical simulations we demonstrate 
that our hypothesis describes the hierarchies of three-frequency resonances in representative dynamical systems. 
We conjecture that this organization may be universal across a large class of three-frequency systems. 
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Nonlinear systems with two competing frequencies show 
resonances or lockings, in which the system locks into a res- 
onant periodic response which has a rational frequency ratio 
[jl|]. The locking increases with nonlinearity, from none in 
the linear regime, to a critical situation where the system is 
everywhere resonant. The subcritical system has quasiperi- 
odic responses between different lockings, while at supercrit- 
ical values of the nonlinearity, chaotic as well as periodic and 
quasiperiodic responses may occur. Resonances have been in- 
vestigated theoretically and experimentally in many nonlinear 
systems, and their distribution in parameter space in the form 
of a devil's staircase is now well understood, from the num- 
ber theoretical concept of Farey trees [§-0] . However, all this 
applies to resonances generated by the interaction of two fre- 
quencies. Far less is known, by comparison, when there are 
three or more interacting frequencies. 

Adding another frequency allows new phenomena to take 
place. Now as well as (two-frequency) resonance as before, 
there is a further possibility: three-frequency resonance, also 
known as weak resonance or partial mode locking. Three- 
frequency resonances are given by the nontrivial solutions of 
the equation afo + bfi + c/2 = 0, where a, b, and c are in- 
tegers, fi and /2 are the forcing frequencies, and fo is the 
resonant response. They form a web in the parameter space of 
the frequencies j^— [TT|] . In this paper, we hypothesize a local 
ordering of three-frequency resonances based on generaliz- 
ing the Farey tree of two-frequency systems to three frequen- 
cies. We perform experiments and numerical simulations to 
show that our hypothesis is justified in representative dynam- 
ical systems with three interacting frequencies: a quasiperi- 
odically forced circle map, a pair of parametrically coupled 
forced nonlinear oscillators, and an experimental system con- 
sisting of an electronic circuit of forced phase-locked loops. 
Our observations lead us to conjecture that the ordering we 
predict may be universal in a large class of dynamical systems 
with three interacting frequencies. 

Firstly, we revise continued fractions and the Farey tree for 
the case of two frequencies. Consider a two-frequency system 
with autonomous frequency fo and external frequency fi . Let 
/ = fi/fo- The aim is to define a sequence of rational num- 
bers that converges to /. Strong convergence [ jl2[ ] is defined 



for rational fractions Pi/qi, (pi,qi) £ Z as 
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for all (pi,qi) for any q { < q n . Given /, p n and q n 
are obtained by expanding / in continued fractions / = 
(ai, 02,03, . . .), and truncating the expansion as p n /q n = 
(ai, 02, 03, . . . a n ) [|l3|]. The p n /q n are then the strong con- 
vergents of /. They give the sequence of fractions with lowest 
monotonically increasing denominators that converges to /. 

The physically motivated hypothesis invoked to explain the 
local ordering of the hierarchy of (two-frequency) resonances 
is that the larger the denominator, the smaller the plateau. The 
fraction with smallest denominator between pjq and r/s, if 
they are sufficiently close that \qr — ps\ — 1, when they are 
called adjacents, is (p+r) / (q+s). This fraction, known as the 
mediant, is then the most important resonance in the interval. 
Repeatedly performing the mediant operation 
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on a pair of adjacent rational numbers, we obtain a Farey tree. 
The Farey tree provides a qualitative local ordering of two- 
frequency resonances, and gives rise to a structure of plateaux 
at all rationals known as the devil's staircase. The devil's stair- 
case in turn is the skeleton for the layout of the resonances in 
parameter space as Arnold tongues . 

Now consider the case of three frequencies, one internal fo, 
and two external f\ and f 2 . We may divide through by the au- 
tonomous frequency f , to give f\ = /1//0, and f\ = f 2 / fo- 
We aim to come up with two convergent sequences of ratio- 
nal numbers with the same denominator, p n /k n and q n /k n , 
which are strong convergents to f\ and f\ respectively. As 
before, strong convergence is defined as 
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Thus (pn/k nj q n /k n ) are best rational approximants if 
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for all triplets of integers (p^, qi,ki) for any ki < k n . So we 
may write 
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which by extension from the two-frequency case we hypothe- 
size to be the largest plateau between fi/p n and f2/q n - Start- 
ing with (/1 , /|) instead of (fi,f 2 ), we obtain the general- 
ized mediant /* = (f 2 - fi) / (q n ~ Pn)\ both f s and /* are 
shown in Fig. [jja). The generalized mediant operation 



(14) 



where we wish to obtain the integers p n , q n and k n . This gen- 
eral problem has not been solved Jl4| , p^[ ], however, we may 
set e\ = e 2 , so that both approximations should be equally 
good or bad. If we do this, we can equate \k n f\ — p n \ — 
\knf2 ~ Qn\, and ask what is k n . There are two solutions 
k n = (Qn i Pn)/(f 2 =^ fl )• At m is point we must remem- 
ber that k n is an integer, so these solutions require that the 
frequencies be rescaled by f\ ± f\ . For which we define for 
the first solution A = f\/{f\ + f\), j 2 = /*/(/* + fl), 
and similarly for the other solution /* = fl/(f 2 — fl), 

/I = fiKfi 
ent e's 



fl). The two solutions give rise to differ- 



£ = \(Pn + q n )fi ~Pn\ = \(Pn + qn)f2 q n \, (8) 

£* = | fen ~Pn)fi ~Pn\ = | (<2n ~ Pn)U Qui (9) 

from which one can obtain e/e* = \ (f 2 — fi)/ (fi + f 2 ) \ < 1- 
So in this sense the (fi , f 2 ) solution is better than the (/* , f 2 ) 
solution. 

Now sticking with the (/i, f 2 ) solution, p n , (p n -\-q n ), and 
q n are obtained from the continued fraction expansions of f\ 
and f 2 . Since 
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if we have the continued fraction expansion of /1//2 = 
(ai,a 2 ,a 3 , . . .), that of fi = {a\ + l,a 2 ,a 3 , . . .), and f 2 = 
(1, ai, <22, a3, . . .) Hence if p n /q n is the nth strong conver- 
gent of /1//2, or equivalently of /1//2, given by this contin- 
ued fraction expansion, p n / (p n +q n ) and g n / (p n + q n ) are 
the strong convergents of fi and /2 respectively. If p n /q n is 
a such a convergent of /1//2, we may define as generalized 
adjacents any pair of (fi/ri, fj/rj), with / £ M and r G Z, 
that satisfy 



= l/l^n ~ f2Pnl 
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From this definition, the subharmonics fi/p n and /2/^n are 
generalized adjacents, and the mediant between them is 
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then provides us with a generalized Farey tree for three fre- 
quencies; m Fig. 0(b) we show the first three levels of the tree 
obtained by recursive application of Eq. ( Jl 4| > to the adjacents 

/i/p n and/ s *. 
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FIG. 1. (a) Schematic diagram of the frequency line, showing the 
relative position of f s and /* and their respective p n , q n harmonics 
that approximate f\ and fi at equal distance, (b) Generalized Farey 
tree starting from the adjacents fi/p n and /*. The first mediant is 
; at the second level we obtain f s and (2/2 — fi )/ (2g n — Pn), 
and so on. 



Let us take as an example a three-frequency system with 
the two external frequencies set to f\ = 2100 Hz and f 2 = 
3600 Hz. The frequency ratio /1//2 is then 7/12. The con- 
tinued fraction expansion for /1//2 is (1,1,2,1,1), and the 
different truncations of this produce the convergents of 7/12, 
which are 1/1, 1/2, 3/5, and 4/7. So we take 4/7 as an approx- 
imation to the higher-order rational 7/12, or equivalently in 
terms of the original frequencies, /1//2 =2000 Hz/3500 Hz 
approximates /1//2 =2100 Hz/3600 Hz with a detuning of 
100 Hz. Between the adjacents fi/4 = 525 Hz and f 2 /7 = 
514.3 Hz lies the mediant f s = (f x + /2VH = 518.2 Hz, 
which we have hypothesized to be the widest resonance in 
this interval, and recursively applying the mediant operation 
gives us the Farey tree that predicts the entire hierachy of res- 
onances in the interval. We now proceed to test this hypothesis 
in three-frequency dynamical systems. 
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FIG. 2. (a) The electronic circuit in block-diagram form. Two 
coupled voltage-controlled oscillators (VCOs) are forced with two 
independent external forces, (b) Three-frequency devil's staircase 
for the circuit. The straight lines correspond to intervals of the DC 
offset within which the fundamental frequency of the output remains 
constant. We varied the DC offset of input 1 in steps of 1 mV in the 
interval 0.90 V to 1.05 V, which corresponds to frequency responses 
of the circuit between 514.28 Hz and 525.00 Hz. Only intervals with 
a stability width greater than 2 mV are plotted. The external frequen- 
cies are fixed at 2100 Hz and 3600 Hz. (c) The generalized Farey tree 
predicts the organization of all the frequency values observed in (b). 

We have constructed the three-frequency electronic oscil- 
lator shown in Fig. ||(a) [jl6|]. Two phase-locked loop (PLL) 
circuits made up of voltage-controlled oscillators (VCOs) are 
coupled through a lowpass and integrator network and forced 
with two independent periodic forces. Typically, the outputs 
of the two oscillators are synchronized 1/1. We filter the out- 
puts of the two oscillators in order to attenuate the compo- 
nents at the external frequencies and we measure the output 
frequencies of both oscillators. If the two values coincide we 



plot them against the DC offset of oscillator 1 , which is used 
as the control parameter (for the circuit a variation in the mean 
value of the ith external force is equivalent to a linear change 
in the natural frequency of the ith oscillator). The results are 
presented in Fig. [z](b) for the interval (/2/7, /i/4) for fi and 
/2 fixed at 2100 Hz and 3600 Hz respectively. This appears to 
be a typical devil's staircase familiar from periodically driven 
oscillators. It is, however, a three-frequency devil's staircase: 
the plateaux correspond to solutions with three linearly de- 
pendent basic frequencies (the frequency plotted plus the two 
forcing frequencies) instead of to periodic solutions. As pre- 
dicted, the generalized mediant f s = 518.2 Hz is the largest 
resonance in the interval, moreover, the generalized Farey tree 
shown in Fig. ||(c) gives the entire hierarchy of resonances of 
Fig. 1(b). 
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FIG. 3. The three-frequency devil's staircase for the system of 
ordinary differential equations of Eq. (|l5|). The system parameters 
are equivalent to those of Fig. ^(b). 

We have integrated an exactly soluble system of ordinary 
differential equations with three interacting frequencies | 

Xi + (4bx 2 - 2a)Xi + b 2 x\ - 2abxf 



+ ((J 0i (t) 2 + a 2 )xi = fi(t) 
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for i = 1,2. The external forces and the coupling term 
are chosen in such a way as to preserve the piecewise in- 
tegrability of the overall system. The oscillators are cou- 
pled parametrically, their intrinsic frequency changing every 
time the coordinate of one of them changes sign cooi(t) = 
cjoi + sgnui(t) sgn Uj(t)Ai. The ith oscillator is driven by 
an impulsive external force fi(t) — Vei J2 n — nr Ei) of 
frequency uoi whose effect is to produce a discontinuity Ve% 
in the oscillator velocity at times tltei = n/ (270^). We have 
performed a power spectrum analysis of the output of both os- 
cillators while varying the intrinsic frequency of oscillator 1 . 
We display in Fig. [3] the most prominent peak in each spec- 
trum against o;oi for a parameter region equivalent to that of 
Fig. ^(b), and in which the two oscillators are also synchro- 
nized 1/1. All detectable resonances are again well described 
by the Farey tree structure of Fig. ||(c). 

Our final example of a three-frequency system is the 
quasiperiodically forced circle map 
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The quasiperiodic sequence ft n is the time interval between 
successive pulses of a sequence composed of the superposi- 
tion of two periodic subsequences, one of period equal to one 
and the other of period Ti /T 2 < 1 (with no loss of general- 
ity), multiplied by the value of the intrinsic frequency of the 
oscillator. In Fig. ||(a) we demonstrate that, for the same input 
frequencies as in the previous two cases, the output in the form 
of a devil's staircase is qualitatively unchanged; once again its 
organization is given by the generalized Farey tree. Moreover, 
Fig. ^(b) shows how three-frequency resonances are arranged 
globally in the form of devil's ramps in the parameter space. 

By generalizing the known Farey tree structure of two- 
frequency resonances for three frequencies we have predicted 
the organization of three-frequency resonances in dynami- 
cal systems with three interacting frequencies. Our results 
for three different three-frequency dynamical systems for the 
same ratio of forcing frequencies fi / show that in each ex- 
ample the predicted generalized Farey tree hierarchy is ob- 
served. We have repeated these experiments and simula- 
tions for different frequency ratios, both rational, as the ex- 
ample presented here, and irrational, for example /1//2 = 
(1 + V / 5)/2, the golden ratio. In every case we have exam- 
ined, the local ordering of resonances around the convergents 
in the devil's staircase — slices through the devil's ramps of 
Fig. @(b) — is well described by the generalized Farey tree hi- 
erarchy. We conclude that the organization we have described 
here is widespread, and conjecture that it may be universal in 
a large class of three-frequency systems. 
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FIG. 4. Quasiperiodically forced circle 

map: (a) Three-frequency devil's staircase for forcing frequencies 
fi = 1, fi — 12/7. (b) Devil's ramps: the rotation number as a 
function of the external frequency ratio and the intrinsic frequency 
shows the global organization of three-frequency resonances. 
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